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Aperiodic tilings are non-periodic tilings characterized by local constraints. They play a key role 
in the proof of the undecidability of the domino problem (1964) and naturally model quasicrystals 
(discovered in 1982). A central question is to characterize, among a class of non-periodic tilings, the 
aperiodic ones. In this paper, we answer this question for the well-studied class of non-periodic tilings 
obtained by digitizing irrational vector spaces. Namely, we prove that such tilings are aperiodic if 
and only if the digitized vector spaces are computable. 

1 Introduction 

A tiling is a covering of a given space by interior-disjoint compacts called tiles. The shape of tiles yields 
constraints on the way tiles can locally be arranged - one speaks about local rules (think, e.g., of the 
bumps and dents of a jigsaw puzzle). Additional local rules can also be set by arbitrarily specifying 
how tiles can be neighboor (tiles can also be available in different colors in order to allow them playing 
different roles). 

Tilings have been studied in computer sciences in the early 60's by the logician Hao Wang, who set 
the so-called domino problem ( lfl4l ): can one decide whether a given finite set of tiles can form a tiling 
of the plane (each tile can be used several times)? His student Robert Berger proved the undecidability 
of this problem (HI). The two key ingredients of the proof are, first, the simulation of Turing computa- 
tions by tilings of the plane and, second, the existence of aperiodic tile sets, that are finite tile sets which 
do tile the plane but only in a non-periodic fashion (Berger explicitly described the first ever such tile set). 

The interest in aperiodic tilings {i.e., tilings by aperiodic tile sets) spreaded beyond computer sciences 
two decades later, when new non-periodic crystals soon called quasicrystals were incidentally discovered 
by the chemist Dan Shechtman ( |[TTI ). The connection with aperiodic tilings was indeed quickly done, 
with tiles and local rules respectively modelling atom clusters and finite range energetic interactions. 
The issue that now concerns theoretical physicists is the classification of all the possible quasicrystalline 
structure, in the spirit of the Bravais-Fedorov classification of crystalline structures. 

A promising approach is the one opened by Leonid Levitov in [9 ]. He considered non-periodic planar 
tilings, that are digitizations of irrational vector spaces, and searched algebraic conditions on vector space 
parameters for the existence of local rules. This approach led to numerous results (ll^l5ll6ll7l[^ [T2llT3i ). 
but no complete characterization of aperiodic planar tilings has yet been obtained. 

The aim of this paper is to move a step forward in the above approach by enriching geometric meth- 
ods with calculability, in the spirit of the first works on aperiodic tile sets. Our main result (Theorem [TJ 
below) is that a planar tiling admits local rules if and only if it is a digitization of a vector space whose 
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parameters are computable. 

Our result thus provides a complete characterization - at least if we do not care about parameters as 
the size of tile sets or the precision of digitizations. Of course, these parameters are important w.r.t. the 
quasicrystal modelization. In particular, the tile sets obtained via the Levitov approach are much smaller 
than the huge ones that we get by simulating computations of Turing machines (the planar tilings ob- 
tained in the former case are however restricted to algebraic parameters, versus computable parameters 
in the latter case). The quest for sharper classifications of aperiodic planar tilings thus remains open, 
with our result nevertheless limiting the horizon of any such classification. 

The rest of the paper is organized as follows. In Section [2j we introduce the formalism which 
allows us to state us our main result, Theorem [T] Section [3] shows that one cannot expect to go beyond 
computability with local rules. The following sections are devoted to prove that local rules allow to 
indeed reach the computability barrier. Specifically, Section [4] introduces quasi-Sturmian words, which 
are particular non-periodic words whose letters are indexed by Z. Local rules cannot characterize such 
words, but Section [5] shows that they can characterize the set of two-dimensional words (letters indexed 
by Z 2 ) whose lines are quasi-Sturmian words. The key ingredient is a result simultaneously obtained by 
[3 ] and [4]. We then show in Section|6]how to transfer this result onto planar tilings which are digitization 
of planes in M 3 , and we finally extend this to any planar tiling in Sec. [7] 

2 Formalization 

n — > d tilings 

Let vi , . . . , v„ be pairwise non-colinear vectors of M. d , n > d > 0. A ;i — > d tile is, up to translation, a 
parallelotope generated by d of the v,'s, i.e., the linear combinations with coefficient in [0, 1] of d of the 
v/'s. Then, an — > d tiling is a face-to-face tiling of M. d by n — > d tiles, i.e., a covering of W 1 by n — > d 
tiles which can intersect only on full faces of dimension less than d. 

Lift 

Let e*i,...,e n be the canonical basis of W 1 . Given a n — > d tiling, we first arbitrarily map one of its 
vertex to £ W, then we map each tile generated by v,-, , . . . ,v; d onto the cf -dimensional face of a unit 
hypercube of Z" generated by e\ d , with two tiles adjacent along an edge v,- being mapped onto two 

faces adjacent along an edge This defines, up to the choice of the initial vertex, the lift of the tiling. 
This is a digital J-dimensional manifold in W, whence d and n — d are respectively called the dimension 
and the codimension of the tiling. 

Computable planar tilings 

A n — > d tiling is said to be planar if there are a (i-dimensional vector subspace Fct" and a positive 
integer w such that this tiling can be lifted into the slice V + [0,w)". The space V is called the slope of 
the tiling and the smallest suitable w its width (both are uniquely defined). A planar tiling is said to be 
computable if its slope is computable, i.e., admits a basis of vectors with computable coordinates that is 
to say they can be computed to within any desired precision by a Turing machine. 

Local rules 

A planar n^-d tiling of slope V is said to admit local rules if, when n—^d tiles are available in colors 
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chosen among a given finite set, the way two colored tiles can intersect can be thoroughly specified so 
that the allowed colored tilings form a non-empty set of colored planar n^d tilings of slope V and uni- 
formly bounded width w. In other words, the slope V is characterized by finitely many rules governing 
the way tiles locally match, with colors allowing each tile to play different roles. In J9), local rules are 
said to be strong if w = 1, weak otherwise. 

We are now in a position to state our main result: 

Theorem 1 A planar tiling admits local rules if and only if it is computable. 

3 The Computability barrier 

Consider W endowed with the norm ||v||oo which gives the maximum of the absolute value of the coor- 
dinates of v. Let S = {x G W : 11x1100 = 1}, for two <i-dimensional vector spaces V and W of W, we can 
define the distance 



Since V and W have the same dimension, then the maximum in the expression above is always attained 
by both expressions simultaneously. The set of cf-dimensional vector spaces is compact by this distance. 
Moreover if we know the computable basis of V and W then d(V,W) is also computable. With this for- 
malism one obtains an equivalent definition of computable vector spaces: V is computable if there exists 
a Turing machine such that on the input n E N it gives a rational basis of a vector space W„ such that 



Proposition 1 If a planar tiling admits local rules, then it is computable. 

Proof. Consider a planar n — ^ d tiling of slope V and width w which admits local rules and take P r the 
set of all the diameter r patterns centred on of colored tilings allowed by these local rules (this takes 
exponential but finite time in r). Let 3£ r be the set of <i-dimensional vector spaces which admit a basis 
given by d vectors associated at a border vertex in the lift of a pattern of P r . The set X r is finite, moreover 
there exists W G 3£ r such that d(W, V) <j. Since for sufficiently large r all vector spaces of 3£ r are near 
of V (if not by compacity one obtains one other slope for the n — > d tiling), to obtain an approximation 
of V with an error bound e, we take the first r>^ such that for all W\,W2 £ 3£ r one has d{W\ , W2) < §. 
In this case, for all W G X r one has d(V, W) < £ so all basis of norm one of W are an approximation of 
a basis of V. Thus V is a computable vector space. □ 



4 Quasi-Sturmian words 

Consider the set {0, 1} Z of bi-infinite words over the alphabet {0, 1} endowed with the metric d defined, 
for any u and v, by 




d{V,W n ) < 



n 



We here show the easiest part of Theorem [T] 



d(u,v) := sup I \u(p)u(p + 1) . . .u(q)\o - \v(p)v(p + 1) . ..v(q)\ 

p<i 



136 



Local rules for computable planar tilings 



where |w|o denote the number of occurences of the letter in the finite word w. In other terms, the 
distance between two words is the maximum balance between their finite factors which begin and start 
at the same positions. 

Define the Sturmian word s pa G {0, 1} Z of slope a G [0, 1] and intercept p by 



Sturmian words have been extensively studied (see, e.g., llTOll for a detailed account). In particular, they 
can be seen as 2 — > 1 planar tilings of width 1 (the digitized vector space is here a line of slope a, and 
the tiles are letters and 1). Classic properties of Sturmian words easily yield: 

Proposition 2 Sturmian words with equal slopes are at distance at most one. 

Proof. Two sturmian words u and v with equal slopes are known to have the same finite factors. Any 
two factors of respectively u and v which begin and start at the same positions are thus also factors of u 
only - at different position but with the same number of letters. This yields the bound 



The words at distance at most one from a Sturmian word of slope a are however not all Sturmian. 
We call them quasi-Sturmian (of slope a). They can be seen as 2 — > 1 planar tilings of width 2. The easy 
following proposition will be useful to link Sturmian and quasi-Sturmian words: 

Proposition 3 Two words in {0, 1} Z are at distance at most one if and only if each can be obtained from 
the other by performing letter replacements — > 1 or 1 — > 0, without two consecutive replacements of 
the same type. 

Proof. Let u and v in {0, 1} Z at distance at most one. Performing on u replacements at each position i 
where u{i) ^ v(i) yields v. If two consecutive replacements, say at position p and q, have the same type, 
then the balance between u{p) . . .u{q) and v{p) . . - v{q) is two, hence d(u,v) > 2. The type of replace- 
ments thus necessarily alternates. 

Conversely, assume that v G {0, 1} Z is obtained from u G {0, 1} Z by performing replacements whose type 
alternates. Given p < q, consider the number of replacements between positions p and q: the balance 
between u(p) . . . u(q) and v(p) . . .v(q) is if this number is even, 1 otherwise, hence d(u,v) < 1. □ 

Since the replacements to transform u in v alternate, their sequence can be encoded by w G {0, 1} Z : 
reading 01 (resp. 10) at position i means that a replacement 0—^1 (resp 1—^0) occurs at position i. Such 
a word w is moreover unique, except if u = v in which case both w = Z and w = l z suit. This coding 
will be used in the proof of Prop. [4| Figures [T] and [2] illustrate this. 

5 A Sofic subshift 

The main result of [3l |H is phrased in terms of symbolic dynamics. Let us first recall this formalism, 
which is convenient to keep, and then explain the correspondance with tilings. Given a finite alphabet 
£/, a configuration is a word indexed by Z". Consider a set of finite patterns F, the subshift of forbidden 



Sp,a(n) = 44> (p +na) mod 1 G [0, 1 - a). 




This bound is known to be at most one for Sturmian words (and only them). 



□ 
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Figure 1: Two (factors of) bi-infinite words with emphasized replacements (first two lines). The corre- 
sponding coding (last line). 
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Figure 2: Two (factors of) Sturmian words with the same slope and emphasized replacements (first two 
lines). The corresponding coding (last line). 



patterns F is the set of configurations where no pattern in F appears. A subshift S is said to have, finite 
type if there is a finite set of forbidden patterns. A subshift S is said to be sofic if there is a subshift of 
finite type S' and a map from the alphabet of 5" to the alphabet of S, called factor map, which maps S' 
onto 5. Last, a subshift is said to be effective if its forbidden patterns can be enumerated by a Turing 
machine. In terms of tilings, a configuration can be seen as a tiling of E" by tiles which are colored unit 
hypercubes. Forbidden patterns then correspond to local rules, and the factor map simply correspond 
to a map on colors of tiles. We are now in a position to use the result of [£3j |4). It is (constructively) 
proven that any effective J-dimensional subshift can be obtained as the projective subaction of a (d+ 1)- 
dimensional sofic subshift, that is, the projection onto d given coordinates. 



First, let Z a be the two-dimensional subshift whose configurations are obtained by copying on each 
row a given Sturmian word of slope a: 

Z a = |«€ {0,1} z2 , BpeR, VmeZ, u(m,-)=s a>p y 

When a is computable, we compute the word s a ,o until to obtain n + 1 distinct factors of size n, thus 
the finite patterns of s a ,p are recursively enumerated. It follows that the projective subaction of Z a is 
effective, thus Z a is sofic when a is computable by 01 Si- Say, e.g., that Z a is the projection onto the 
first coordinate of a finite type subshift Z a of ({0, 1} x B), for some finite alphabet B. 

Then, extend Z a by the subshift Z' a whose elements have rows at distance at most one from the 
Sturmian word of slope a and intercept 0: 

Z' a := |mg{0,1} z2 , Vm6Z, d(u(m,-),s afi ) < lj. 

Prop. [2] indeed yields Z a C Z' a and that the choice of the intercept - here - has no importance. Let us 
constructively prove: 

Proposition 4 The subshift Z' a is sofic when a is computable. 

Proof. Let Xi u ...j k denotes the projection on the ii-th,. . . ,4-th coordinates. Let Z' a be the subshift of 
({0, 1}xBx{0,1}) z2 such that u G Z' a if and only if 7ln(u) £ Z a and, for any (m,n) G I?: 

Ti-i{u(m,ri)) < K^(u(m,n + 1)) =£- 7Z\(u{m,n)) = 0, 
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Ki{u{m,n)) > 713(11(111, n + 1)) K\(u(m,n)) = 1. 

Clearly, Z' a has finite type when so does Z a . Now, we claim that the factor map % defined as follows 
maps Z' a onto Z' a : 

, w s f Tl\(u(m,n)) if m(u(m,n)) = Ki{u(m,n + \)) , 
n(u)(m,n) = < ; , . 

[ 1 — Tii(u(m,n)) otherwise. 

Let hGZJ, and fix m G Z. By definition of Z„ and Z a , %\ (u(m, • )) = Sa.p. One thus also has 7l(u(m, •)) = 
s a p , except at each position n such that the two bits 7l-s(u(m,n)) and 7ii,(u(m,n + 1)) differ. At these posi- 
tions, 7t(u(m, ■ )) is obtained by performing on s a<p a replacement of type 7[j,(u(m,n)) — > 7F3 (u(m,n+ 1)). 
The type of these replacements alternate - as the bit runs do - and Prop. [3] yields d(n(u(m, ■ )), J«,p) < 1- 
This shows that 71(G) is in Z' a . Hence, 7l(Z' a ) C Z' a . 

Conversely, let u G Z^. Fix m G Z and choose v G Z a such that %\ (v(m, • )) = s a $- By definition, 
d(u(m,- ),s a ,o) < 1, so we can consider w,„ the coding of the replacements which transform s a .o into 
u(m, ■ ) (see end of Section|4]). Consider u G ({0, 1} xfix{0, l})^ 2 defined by u(m,i) = (v(m,i),w m (i)). 
The way n has been defined yields u G Z^ and 7r(«) = w. Hence, C 7c(Z' a ). □ 

Let us mention that the two-dimensionality plays a fundamental role in the result obtained in 0141, 
hence in the soficity of Z a . It is thus also fundamental in the proof of Prop. |4j although lines seems to 
be there only independantly considered. The interplay between the lines of Z' a is indeed "hidden" in the 
alphabet B and in the forbidden patterns of Z a . 



6 Dimension two and codimension one 

Consider a computable 3—^2 planar tiling defined over vectors vi, V2 and i%. Its slope can be defined by 
its normal vector, say (1, a, j3), which is computable. 

As mentioned above, any sofic subshift can be seen as tilings by a given tile set and local rules. Here, it 
is convenient to model local rules by coloring edges of tiles and assuming that two tiles can be adjacent 
only on edges which have the same color (this is actually the definition used in [14]). Let thus z a and 
Tp be such tile sets, with moreover label or 1 on tiles, such that the tilings by z a and Xp respectively 
correspond, when considering only labels, to Z' a and Z^. 

Let us derive from T a a tile set %' a as follows (see Fig. [5]): 

• each tile in z a with label is sheared along V3 to give in z' a a rhombus tile defined by V2 and V3 
(edge colors are unmodified); 

• each tile in z a with label 1 is sheared along V3 to give in x' a a rhombus tile defined by vi and V3 
(edge colors are unmodified); 

• a rhombus tile defined by vi and V2 is in r' a if and only if its edges have colors which appears on 
vi- or v^-edges of square tiles of z a , with the restriction that two edges meeting at vi or V2 must 
have the same color (transfer tile). 

The idea behind the definition of r' a is simple. Let us call Vi-ribbon of a 3 — ^ 2 tiling a maximal 
sequence of tiles, with two consecutive tiles being adjacent along an edge v,-. Then, one easily sees that 
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Figure 3: The v,'s (left), the shearing of tiles in r a (center) and the transfer tiles (right). 
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Figure 4: A part of an element in Z' a (left) and its embedding into the V3 -ribbons of a 3 — > 2 tiling (right, 
with thick lines drawn on "transfer" tiles indicating the way tile decorations travel between neighbor 
V3 -ribbons.) 



z' a forms the 3 — > 2 tilings whose V3-ribbons embed configurations in Z' a , with the transfer tiles just 
carrying colors between ribbons (see Fig. [4]). 

By proceeding similarly up to a circular permutation on the v,'s, one derives from Zp a tile set 
which forms the 3 — > 2 tilings whose iVribbons embed configurations in Z'n . Let us finally define the tile 
set l' a p as the cartesian product of z' a and , that is, to each pair of identically shaped tiles in x' a x 
corresponds a tile in i' a ^ which has the same shape, with each edge having a color which encodes the 
colors of the pair of corresponding edges in x' a and . The tilings of z' a ^ are thus the 3 — > 2 tilings 
which embed Z' a on their V3-ribbons and on their iVribbons. This allows only planar 3—^2 tilings 
of slope (l,a,j6) and width at most 4 (this can be seen by decomposing any path between two points 
in a V3-ribbon followed by a iVribbon). This moreover allows at least the planar 3—^2 tilings of slope 
(l,0l,j8) and width 1. We thus (constructively) proved: 

Proposition 5 Any computable planar 3—^2 tiling admits local rules. 



7 Higher dimensions and codimensions 

The last step to prove Th. [T] is to extend Prop. [5] to n — > d tilings. Although technical, this last step 
requires no new ideas. 

For higher codimensions, we proceed by induction. Our induction hypothesis is that any effective 
planar n — > 2 tiling admits weak local rules. This holds for n = 3 according to the previous section. Let 
now & be an effective planar (n + 1) — >■ 2 tiling. For any basis vector e,-, we project the lift of & along 
e*i to get the lift of an effective planar n — > 2 tiling, say By assumption, ^ admits local rules: let T; 
be a tile set whose tilings are at distance at most w from 3^. We complete T, by adding the tiles with 
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a v,-edge (that is, the tiles which disappeared from 3? by projecting along with each of these tiles 
having no decoration on its v,-edges, and on the other edges a unique decoration that could be any of 
those appearing on an edge of a tile in T,-. These new tiles thus just transfer decorations between the tiles 
of 3?i (see Fig. [5]for n + 1 = 4). Last, we define the tile set z as the cartesian product of all the T;'s (as we 
did for T a and Zp in the previous section). This allows only n + 1 — > 2 tilings at distance at most W from 
3* - in particular 3T itself. This shows that 3* admits local rules. 




Figure 5: A 4 — > 2 tiling by shaded T, tiles and white additional tiles, with v; being here horizontal. The 
white tiles simply transfer horizontally the decorations of the shaded tiles. By contracting each v,-edges 
to a point, the white ribbons disappear: we get a 3 — > 2 tiling by at distance at most w' from 3^. 



For higher dimensions, we also proceed by induction. Our induction hypothesis is, for a fixed n, that 
any effective planar n—> d tiling, d <n, admits local rules. This holds for d = 2 according to the above 
paragraph. Let now 3* be an effective planar n — > (d + 1) tiling, with d + 1 < n. Fix i G { 1 , . . . , n}. For 
two tiles T and T of ST , write T ~ T' if these tiles share a v,-edge and let ~ be the transitive closure 
of the relation ~. Denote by {3?k)kez the equivalence classes of ~, such that, for any k, 3\ and 5\+\ 
can be connected by a path which does not cross any other equivalence class (the ^'s play the role 
of v,-ribbons in the previous section). By contracting all the v,-edges of a 5\ {flattening), one gets a 
planar n — > d tiling. Its slope moreover depends only on the slope of 3? , and in particular it is effective. 
This allows to see 3T as a sequence of "stacked" parallel effective planar n — )■ d tilings (namely the 
flattened 3\\), with the remaining tiles containing no v ( -edge. By induction, there exists a finite tile set 
%i whose tilings are at bounded distance w from any of the flattened J^'s (since they are all parallel). It 
is straighforward to "unflatten" T, to get a tile set f; whose tilings are at bounded distance w from any of 
the We complete f ; - by adding the tiles without v,-edge (that is, the tiles lying between the stacked 
=5£'s), with decorations being just transferred between consecutive 3%'s along the direction v; (as done 
in the previous section to transfer decorations between consecutive v,-ribbons). The last step is (as in 
the previous section again) to define the cartesian product z of the tile sets % i = 1, . . . ,d: its tilings are 
those at bounded distance w from 3* - in particular 3* itself. This shows that 3* admits local rules. 
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